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ABSTRACT 


The hybrid model of the Reynolds-stress turbulence closure, which was 
developed In the earlier work Of this project, is tested for the computation 
of the flows over a step and disk. Here It Is attempted to Improve the 
redistributive action of the turbulence energy among the Reynolds stresses. 

By evaluating the existing models for the pressure-strain correlation, better 
coefficients are obtained for the prediction of separating shear flows. 
Furthermore, the diffusion rate of the Reynolds stresses Is reevaluated 
adopting several algebraic correlations for the triple-velocity products. The 
models of Cormack et al., Daly-Harlow, Hanjal Ic-Launder, and Shir were tested 
for the reattaching shear flows. It was generally observed that all these 
algebraic models give considerably low values of the triple-velocity 
products. This is attributed to the fact that none of the algebraic models 
can take the convective effect of the triple-velocity products into account In 
the separating shear flows, thus resulting in much lower diffusion rate than 
Reynolds stresses. In order to Improve the evaluation of these quantities 
correction factors are introduced based on the comparison with some 
experimental data. 
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constants used in pressure-strain correlations 

constants used in pressure-strain correlations 

coefficients used in turbulence model 

coefficient tensor 

diameter of the c.vlinder 

diameter of the disk 

diffusion of iijuj 

generation of turbulence kinetic energy 
generation of Reynolds stresses 
step height 

turbulence kinetic energy 
pressure fMictuation 
mean pressure 

generation rate used in pressure-strain rorrelati 

radial coordinate 

raoius of the disk 

fluctuating velocity 

mean velocity 

mean velocity at the inlet 

free stream velocity 

fluctuating velocity in y-direction 

mtan velocity in y-directicn 


X Cartesian coordinate 

y Cartesian ccordinate 

Y height of inlet flow section 

0 

Q, coefficient used in Eq. (22) 
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energy dissipation rate 
dynamic molecular viscosity 
density 
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1. INTRODUCTION 

This report summarizes the study of the behavior of the second- and 
third-moments In the wake region behind bodies. The flow In such a situation 
has been Investigated by formulating the transport equations for the kinematic 
Reynolds stresses iTju'jJ In evaluating the Reynolds stresses in the 
shear region It Is Important to understand the mechanism of energy transfer 
which occurs among the Reynolds-stress components. Separated flow generates 
turbulence resulting In an Increase In the heat and mass transfer rates which 
could be advantageous In some cases and, therefore, preferable (e.g.. Internal 
combustion engines). On the other hand, separation causes pressure losses 
which are undesirable but are sometimes difficult to avoid In aeronautics. 
Whether It is required to create, control or avoid separation and recircula- 
tion, a better understanding of the turbulence phenomena In the situation will 
provide a good basis for attempts to Improve the performance of fluid machin- 
ery which operates creating such separating turbulent flows. 

The recent development of computer techniques for solving the partial 
differential equations governing these flows has made theoretical studies 
easier and less costly by eliminating the need for extensive experimental work 
during design stages. 

A model of the Reynolds stresses was first proposed by Rotta", and has 

3“*8 

been developed and Improved by a number of researchers. The correlation 
of pressure-strain was proposed by Naot et al.^"^ and Launder ct a1.® 

Naot et a1. evaluated the pressure-strain correlation term by integrating over 
space after Inserting two-po1nt correlation functions, whereas Launder et al. 
obtained the results by assuming a fourth-order ten’.or consisting of linear 
Reynolds-stress elements. As a simplified model of the Reynolds-stress 


lil 


closure, an algebraic stress model which does not possess both convection and 

9 

diffusion terms was developed In a similar manner by Rodl. 

In the authors' previous paper^® several proposed closures for the 

Reynolds-stress model were tested for heat transfer characteristics along the 

walls of axisymmetric sudden expansion pipes. After the computation with 

several different models, it was found that the pressure-strain correlation 

0 

proposed by Launder et al. showed slightly better agreement with 
experimental data for Nusselt number distributions along the pipe walls than 
with other models. Moreover, Incorporation of the wall correction terms in 
the pressure-strain also improved predictions by about 5 to 10 percent. 

Here computations of the turbulent flow are made in the region beyond a 
disk by employing the second-order closure model of turbulence. The results 
are compared with the existing experimental data at several different 
steanrwise locations. 

Particular attention Is paid to the evaluation of the pressure-strain 
correlation which plays an important role in transfering turbulence energy 
among the Reynolds-stress components In the separated shear flow region. 

First, the models of Naot et al.^”^ and Launder et al.® are examined, and 
the results computed using these models are compared with existing 
experimental data. Moreover, a simpler model is formulated that provides 
better profiles of the Reynolds stresses in the wake region. 

The study Is further extended to the evaluation of the diffusion action 
of the Reynolds stress. As was discussed in the previous report, the 
change in the triple-velocity products is significantly large in the wake 
region resulting in a considerable variation in the diffusion rate of the 
Reynolds stresses. Thus, it is important to reevaluate the existing models of 
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the third-order closure for better understanding of such diffusion processes 

In the separated shear layers. In the present report four models of the 

third-order closure are examined, and the results are compared with the 

12 

experimental data of Chandrsuda and Bradshaw. The models considered are 

3 13 14 

those proposed by Daly and Harlow,'’ Hanjallc and Launder, Shir,” and 
Cormack et al.^^ Mathematical formulations, along with the results and 
discussion, are presented In subsequent sections. 


2. MATHEMATICAL TORMULATION 


The steady, two-dimensional form of continuity and momentum equations are 
given as: 

Continuity Equation: 


(rUj) ^ 0 


Homentem Equation: 


3 3P a aUi 

iTj ^ 

2.1 Reynolds-Stress Equations 

The transport equations for the Reynolds stresses are given as 


iT (ViV “ <^1j - "ij ^ <‘ij ^j.w “lj 


ali^ an, 

5;^ + “i"k ii;> 
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= (76 - 8 b^)/105 
•= (22 + 64 b^)/105 
a, - 2 (9 + 72 b,)/315 

^ (1! 
= (bg + 8)/11 

aj «• (0 b^ -2)/11 

a^ - (30 -2)/55 


The first model was developed by removing the Isotropic constraint from 
the double velocity two-point correlation tensor. The second model, an exten- 
sion of the first, was proposed in terms of the single point tensor and three 
arbitrary functions of the absolute value of the two-point separation distance 
In this model the value of b^ was determined as 



(!' 


where 

K(r) = €2211 - C ^^22 

and where C,.^ represents the coefficient tensor that Is used to 

J K Xlil 

correlate between the single-point double-velocity and the two-point double- 
velocity products such as 


U,(x) u,(x+r) = (x,r) u (x) u (x) 


(1 
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Tiie third model was obtained by approximating the pressure-strain corre- 
lation by an arbitrary fourth-order tensor where the single-point double- 
velocity correlations v;ere used. 

2 . 3 Third-Order Closure Model 

As is discussed in the last part of the Introduction section, the triple 


velocity fluctuation products, u^UjUj^, should be evaluated appropriately 
in a separating and reattaching flow region. Here the four existing models of 
the third-order closiife, all in algebraic forms, are reviewed and tested for 
the reattaching shear layers. These models are summarized in this section. 

3 

Daly and Harlow obtained an algebraic expression with a simplificaticn 


of such a transport equation of u^UjUj^ as follows 




du.u^ 

1 ^ 


-2-OSrVi ax 


J3 


(19) 


Hanjalic and Launder obtained the following algebraic equation. 




auj^u^ 


au_^. 


Shir^^ proposed the following expression by approximating the first 
Reynolds stress that appears in fq. (19) to be isotropic. 


( 20 ) 


“l“j“k 


- 0.0/1 C„ ^ --rj 

r c oX, 


? 3U.U, 


(21) 


Connack et al. obtained algebraic expressions by approximating the 
experimenlally determined profiles for d.Uj, c and 3 






polynomial with coefficients chosen to give a least-squar*? fit to the data and 
discretized the cross-stream coordinate In earh flow by using the most general 
model for the triple velocity correlation teisor as generated using the ns, ap- 
totic approach of Lumley and Khajeh Nouri.'^ The model they obtained Is as 


follows: 




(^ik 


* *ij^kl 


aP ®^1t. 

* SC * “4 


where 


" "i“j ■ 3 


Out of the twenty parameters, q-, that they had started out with, 
they were able to determWie and optimize the most significant four coeffi- 
cients. The values of these parameters have been recormended for various 
kinds of flows along with the universal value applicable to most of th» flows. 


e.l« X 10 


- 1,72 X 10~'- 


= - 4.00 X lU ' 


0.102 
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In the original forms of (19)-(22), the coefficient Cp Is unity, 
b"t the best values for Cp were Investigated for tlie reattaching shear flows 
by comparison with experimental data. The recommended values are given In 
Sec. 9.2. 


3. NUMERICAL MODEL 


3 . 1 Nume rical Method 

The solution method of the transport equations described In the preceding 
section Is the same finite volume method as the one used In [1] In which the 
differencing scheme Is the modified hybrid scheme of Amano.^^ This scheme 
has a combined mode of convection and diffusion which Is derived by expanding 
the analytical one-dimensional solution up to the fourth-order term. The cell 
structure for mean-velocity components Is the staggered system In which the 
locations of the mean velocities U and V are a half-cell shifted In x- and y- 
directions, respectively. All the normal Reynolds stresses (u^) are 
evaluated at the scalar node point along with P, k and c. However, the 
shear Reynolds stress (uv) Is located at the southwest corner ">f the scalar 
cell. This Is because the main driving strains for the shear stress are 
aU/3y and aV/ax which can easily be evaluated without any Inter- 
positions (Fig. 1). 

3 . 2 B oundary Conditions 

For the examination of the pressure-strain correlation, n cylindrical 

coordinate vjjs used since all the components of the Reynolds •tresses (u . 

~2 2 

V , w , and uv) for two-dimensional problems can be Investigated Instead 

~2 ~~2 

of just two components of normal stresses (u and v ) and the slicar com- 
ponent (uv). The flew domain is shown In Fig. 2. For the Investigation of 
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the triple-velocity products, however, the plane two-dimensional coordinate 
was employed due to the limitation of the experimental data for *hese quan- 
tities. (see Fig. 3) 

There are three different types of boundary conditions to be specified 
for the computation of the flows in Figs. 2 and 3: Inlet, outlet, and wall 

boundary conditions. At the Inlet all fte quantities are specified according 
to the fully developed rondltlon. At the outlet a continuative flow condition 
is applied where gradients of flow properties In the flow direction are zero 
(Neumann conditions), i.e., 34«/3x = 0, where «t> = U^, k, c, u^Uj, etc. 

This outlet is located about 80-120H downstream from the step so that Its 
Influence on the main flow region is negligibly small. 

At the wall boundaries the velocities and turbulence quantities must be 
specified functionally according to the law of the wall. The velocity com- 
ponent normal to the wall is simply set as zero. The wall boundary values for 
k are determined by means of wall functions based on the assumption of a 
logarithmic near-wall velocity distribution w: ith allows the wall-shear stress 
to be extracted from the *,og law* and the vaiue of velocity parallel to the 
wall to be computed along the grid line closest to the wall. The enc'gy dis- 
sipation rate, c. Is evaluated under local equilibrium condition- Near-wall 
effects on the turbulence structure, associated with steep velocity varia- 
tions, are also taken Into account by introducing appropriate modifications to 
the generation and dissipctlon of the turbulence energy and the energy dissi- 
pation rate for the finite volume adjacent to the wall.^^ 

The boundary values »or the Reynolds stresses are determined as 
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- 1.21k 
•= 0.24k 

1 «= 0.55k 
w 


(25) 


-uv = -0.24k + ^ ^ 

p QX 

1n the wall adjacent numerical cells. The details of the derivation of (25) 
is given in ref. [11]. 

4. RESULTS AND DISCUSSION 


4 • 1 Flow Beyond a Disk (Evaluation of Pressure-Strain Correlations) 

The flow field considered here is shown in Fiq. 2. The diameter of the 
cylinder is 1/3 of the disk diameter. The computational domain is subdivided 
into many cortrol volumes in such a way that the sizes increase at the rate of 
2 ~ 3 X in both x and r directions. In this way the size of the numerical 
control volume is reasonably fine in the recirculating region but is rela- 
tively coarse in both the outer region and the downstream region where flows 


are mainly parabolic. 

Several cases with different mesh sizes were tested and the size of 42 x 


42 demonstrated a grid independent state for the flow domain consisting of 
four radii of the disk in height and twenty radii in length. 

Figure 4 demonstrates the axial velocity profiles at several different 
locations behind the disk. Agreement between computed results and the experi- 
mental data of Smyth^^ is very good (within 15X discrepancy). It was also 
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discerned that the velocity profiles were almost Independent of the pressure- 
strain correlation although the computed velocity profiles using all the 
models are not shown In the figure. 

~2 ~2 

Figure 5 shows the nonnal stress profiles of u and v at four dif- 
ferent locations downstream of the disk. The computations v«re made by using 
Models 1, 2 and 3 as defined In Sec. P.2, and the results were compared with 
the experimental data of Smyth. When the coefficient C. •= 2/3 Is used 
for Model 1 (which was recommended by the originators^), the results are all 
similar to each other. This Is mainly because, although all the models are 
derived o.fferently, their final forms are very close. In particular. Models 
2 and 3 113'“ almost the sane values for the coefficient of the first terms on 
the right-hand side of both Eqs. (10) and (11). Naot et al.** obtained the 

o 

value of -0.5 for the coefficient b^ and Launder ct al. chose the value 

0.4 for which results In a^ « a^ « 0.76. Also, It Is noticed 

that the second and the third terms of Models 2 and 3 have a minor Influence 

on the levels of turbulence stresses; thus, the pressure-strain correlation 

may be represented by the first term alone. Note that Model 1 has this term 

only. 

In Fig. 5 the computed results using Model 1 with C 0.4 are also 

♦ 

Q 

shown. This value was recommended by Launder et al. In reference 1o equi- 
librium shear flows. The computed results using the coefficient = 0.4 
show that the levels of u Increase about 20% while those of v decrease 
about the same amount; this Indicates that the redistributive ..'Ction Is 

reduced by about 20;t by changing C. from 0.667 to 0.4. Therefore, the 

4 > 

2 

turbulence energy created due to mean strain of the main stream flow (u ) is 

~~2 

not transferred completely to the normal components (v and w“) of the 
turbulence stresses. 


r 
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It was generally observed that all the models underpredict the levels of 
the Reynolds stresses by 10-50X. In order to predict the Reynolds stresses 
more accurately In recirculating flows, the original models tested for equili- 
brium shear flows need to be revised. For ease of analysis Model 1 was chosen 
to refine prediction of the turbulence level. Two computations were performed 
by using values of 0.2 and 1.2 for C^. The results are shown In Fig. 6. 

As Is depicted In this figure, the small value of coefficient results 

2 2 

In appropriate levels of u but unacceptably low levels of v . On the 

other hand, the large value of gives reversed results; that Is, the 

2 

levels of V are substantially high bsic those of u are too low. 

When the flow reattaches on the center cylinder, the flow starts acceler- 
ating In the downstream direction which causes a high mean nom’al strain as 

~2 

the flow recovers. Thus, the term with u aU/ax In the pressure-strain 
~~2 

correlation In the u -equation becomes higher than those in the fully 
developed flov/s which results In significant energy transfer from u to 
V . This effect must be suppressed to some extent In the present flow 

situation, vjhereas the corresponding component v;hlch balances the energy level 

2 2 
of V is relatively small. Therefore, the redistributive action in the v 

equation needs to be promoted more. 

For the reasons discussed above, the case v/hen the coefficient C. for 

<p 

2 2 
u was decreased to 0.2 and when that for v was raised to 1.2 Is demon- 
strated In the same figure. The results obtained with this treatment show a 
large Improvement (sec Fig. 5). This observation suggests that the 
coefficients of the isotropic generation rates for the Reynolds stresses (l.e. 
the first terms of Fqs. (9) - (11)) should be adjusted by their values by the 
strength of the meen strains, since the flew patterns are strongly affected by 
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the frean strains. The strain variations are particularly complex in reattach- 
ing shear layers being accompanied by recirculating flows. In this way the 
Reynolds stresses can be more accurately evaluated by using such simpler 
formulations. 

In Fig. 7 all the Reynolds stresses computed using the method mentioned 
above are displayed at several different axial locations. 

4.2 Flow Beyond a Step fFvaluation of the Triple-Velocity Correlations) 

Before computing the triple-velocity products in the region behind a 
step, each component of the Reynolds stress equation (3) was solved along with 
the momentum, turbulence energy, and turbulence energy dissipation rate 
equations. 

To establish an optimum grid system, exploratory grid tests were 
performed using 32x32, 42x42 and 52x52 mesh sizes and several different grid 
expanding factors in the downstream direction. Because computations with 
coarse mes i sizes fail to provide a decipherable trend ir the Reynolds stress 
distributions and because the triple-velocity products vary steeply in the 
shear layer region the mesh size of 52x52 was used with an axial length of 50 
step heights for a step ratio of Y^/H = 2.5. The grid expansion factors 
used in the streanrwise and transverse directions were 1.01 and 1.02, 
respectively. 

Figure 8 shows the v'elocity profiles at two different locations behind 
the step. The results aie compared with the experimental data of Chandrsuda 
and Bradshaw^^ obtained with two measurement techniques: a hot-wire, and a 

pressure-probe. Although slight disagreement appears near the wall at x/H = 
5.4, the computed results are within 20 % of the experimental data. 


2 2 

Figure 9 shows the normal Reynolds stresses (u and v ) and the shear 
Reynolds stress (uv) distributions at three different locations In the flow 
field beyond a step. As Is shown In this figure all the Reynolds stresses 
decay In the region downstream of the reattachtnent. In general, agreement 
between the computed results and the experimental data Is fairly good. Based 
upon this agreement with the experimental data, the values of these Reynolds 
stresses were used to evaluate the algebraic models of the triple-velocity 
correlations (Eqs. (19) - (22)). 

Figure 10 represents the distribution of the triple-velocity products at 
four different locations In the reattaching shear layer. It is commonly 
observed that all the models underpredict the levels of triple-velocity 
products everywhere. While the model of Shir consistently gives low levels, 
the model of Cormack et al. gives the highest levels for uvv and Tvv. The 
model of Oaly-Harlow, which is most frequently used gives relatively lower 
levels e.'^cepl for uTiv, whereas the model of Hanjallc-Launder gives reasonably 
high levels for every component of the triple-velocity products. 

In order to adjust for the discrepancy In the peak va/ues of the 
triple-velocity products, the multiplying factor Cp In Eqs. (19) - (22) was 
Incorporated to suitably compensate for the difference In the levels. This 
was done by dividing the experimental peak values of u^Ujuj^ by those 
predicted by the particular algebraic model. These factors were obtained for 
uuv, uvv, and Tvv profiles at x/H locations of 5.4, 8.4, 10.3, and 12.3 for 
all four algebraic models. 

lable 1 shows the values of these factors for Individial components and 
models. The overall factor for a particular model is then obtained by 
averaging all the values of Cp for all the components or ii^UjU^ for 
that model. 
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Table 1. Correction Factors for u^UjUj, for 
Different Algebraic Models 


Models 

uuv 

vvv 

uvv 

Overall 

Average 

Oaly-Harlow 

2.66 

4.93 

6.45 

4.68 

Hanjalic-Launder 

4.18 

3.74 

4.76 

4.22 

Shir 

7.68 

14.20 

18.66 

13.51 

Cormack et al. 

- 

2.12 

3.52 

2.82 


The triple-velocity products were recomputed by using the overall 
averaged factors and the results are shown in Fig, 11. It is interesting to 
note that the peak values obtained by using the model of Hanjalic-Launder 
agree with the experimental data for all the components, whereas those 
obtained by other models agree with the experimental data only for uvv 
component. This indicates *’iat, although the l-'vels of the triple-velocity 
products can easily be adjusted according to flow conditions, we cannot 
necessarily obtain universally improved results for every component of 
li^UjUj^ when using models other than th.at of Hanjalic-Launder. 

The reason that the model of Hanjalic-Launder gives better results for 
all the components of tT^UjFj, Is because in the derivation process of 
this model the generation rates doe to turbulence stresses as well as the 
diffusion rates of u^Uj ware taken into account, even though the 
convection and the generation rates due to mean strains were neglected (see 
ref. 13). The generation caused by the stresses is particularly significant 
in the reattaching shear layer. 



Unlike this model the models of Oaly-Harlow and Shir do not include 

non-homogeneous flow effects. That is, when the Shir model, for example, is 

~2 

used, the transverse gradient of u is the only component that evaluates the 

behavior of the triple-velocity product, uuv, in the transport equation of 

u . The model of Oaly-Harlow has an additional component that accounts for 
2 

the variation of u in the streamwise direction for the corresponding 

product. However, these two models are notably limited in comparison with the 

Hanjalic-Launder model which accounts for variations of the shear stress in 

~2 

both directions in addition to the u variations. It should, however, be 
noted that all the models would produce similar results in an Isotropic 
turbulence flow field. 

The model of Cormack et al. was evaluated only for vvv and uvv since the 
values computed for uuv using its original form were not high enough to 
estimate the factor Cp within a reasonable range (see Fig. 10). 

Finally, the triple-velocity products were recomputed by employing the 
individual correction factor Cp for each component and the results are shown 
in Fig. 12. in this case all the models agree with each other due to their 
component-wise correction factors. 

5. SUMHA1IZ1NG RFMARKS 

Pressure-Strain Correlation s 

6 8 

1. The models by Naot et al and Laundc" et al. effectively give 
similar results and both models are reasonably reliable. 

2. The energy redistribution can be improved for reattaching shear 
flows by taking the effects of moan strain into account. 



I 



Triplc-VcIocUv Correlations 

1. The triple-velocity producis need to be predicted accurately in the 
reattaching siicar layer in order to evaluate the diffusion process of 
the Reynolds stress appropriately. The behavior of the 
triple-velocity products in such a complex turbulent flovj is 
different from that in sittpler flows. 

2. All the existing algebraic models for the triple-velocity products 
underpredict the levels of u^UjUj^ in the reattaching shear 
layer. The predicted levels can easily be improved by using a 
correction factor. 

3. With the exception of the model of Hanjal ic-Launder, all models 
cannot Improve the prediction simply by employing a single value for 
Cp. Thus, it is difficult to obtain a unique value for the 
correction factor. This is because the Hanjaiic-Launder model is the 
only one that includes the generation terms due to Reynolds stresses. 

Finally, it was observed that none of the above models accurately 
predicts the overall levels of the triple-velocity product-.. This is 
primarily attributed to the fact that the convection effect of the 
triple-velocity products is never taken into account. This effect may be 
small in simple shear layers but is significant when the shear layer 
reattaches on a solid wall transfering high convective rates of u^UjUj^’ 
into turbulence diffusive energy near the wall and, subsequ»-ntly, these high 
levels are transported in the downstream direction. In con*, equenre, a model 
that can take into account this process must be developed end tested for 
complex turbulent flows. 
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Fig. 3. Coriiputational domain for the step Hew. 
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Fig. 5. Reynolds-stress profiles behind the disk 
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Fig. 12. Triple— velocity product profiles beyond the step (Individual 
Lorrectlon factors ore used for each component). 
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